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Resumen 


Considerando que el esfuerzo cortante total en el flujo del agua resulta de la suma 
de un esfuerzo viscoso y uno turbulento, se introduce una longitud de mezcla 
hiperbólica en la hipótesis de Newton para describir la relación entre el esfuerzo 
cortante total y la derivada del perfil medio de las velocidades. La longitud 
hiperbólica se reduce a la longitud de mezcla de Prandtl, que es proporcional a la 
distancia a la pared cuando ésta es grande. La existencia de un punto de inflexión 
en el perfil de las velocidades induce un valor teórico de la constante de von 
Kármán muy cercano a su valor experimental. Con este perfil se deducen fórmulas 
para las velocidades medias del fluido en un tubo y en un canal rectangular; de 
éstas se deducen las leyes de resistencia al flujo. La ley de resistencia en el tubo 
llevada a la forma de Darcy-Weisbach indica que el factor de fricción presenta los 
comportamientos extremos de Poiseuille y de Prandtl. La rugosidad de la pared a 
altos números de Reynolds se introduce siguiendo la idea de Colebrook y White, 
al adicionar a la escala de la turbulencia otra proporcional al tamaño de los granos 
de arena de Nikuradse; el factor de fricción tiene el mismo comportamiento que el 
mostrado en el diagrama de Moody. La escala dependerá del número de Reynolds 
en régimen transitorio de acuerdo con los resultados experimentales de Nikuradse, 
y su dependencia a pequeños y moderados valores del número de Reynolds y de la 
rugosidad es un sujeto de estudio. 


Palabras clave: un valor teórico de la constante de von Kármán, efectos de la 
rugosidad de la pared. 


Introducción 


En hidráulica de canales y tuberías, hidrología 
superficial, ingeniería de riego y drenaje, entre 
otras disciplinas, es de fundamental importancia 
la utilización de una ley de resistencia al flujo 
para estimar la fuerza de fricción. La ley de 
resistencia relaciona la pendiente de fricción (J) 
con la velocidad media (U) y el radio hidráulico 
(R,) (Du Buat, 1786; Trueba-Coronel, 1954). 

En flujo a superficie libre es común utilizar 
la ley en potencia U = cR',J”, que incluye las 
leyes extremas en flujo turbulento de Chézy 
(1776) (b = 1/2, d = 1/2) y en flujo laminar de 
Poiseuille (1840, 1841) (b=2, d=1) originalmente 
obtenida en tubos capilares. En un estudio 


sobre el acoplamiento de las ecuaciones de 
Saint-Venant (1871) y de Richards (1931) en el 
riego por gravedad en un canal rectangular 
muy ancho, donde R, = h, siendo h el tirante 
de agua, Fuentes et al. (2004) establecen que 
la singularidad presente en el tiempo inicial 
requiere que los exponentes estén relacionados 
por b = 3d — 1. Esta relación es exactamente 
satisfecha por las leyes de Chézy y Poiseuille, 
aproximadamente por la ley de Hazen-Williams 
(b=0.63, d =0.54) utilizada en el diseño de redes 
de tuberías (King et al., 1952) y claramente la ley 
de Mamning (1891) (b=2/3, d =1/2) es excluida 
por la singularidad. 

La ley de resistencia en potencia que 
satisface la singularidad, haciendo intervenir 
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el coeficiente de viscosidad cinemática (v) y la 
aceleración gravitacional (<), puede ser escrita 
como U = kv'4p3+(g)) donde 1/2<d<1yk 
es un coeficiente empírico; la expresión para 
el caudal unitario es q = Uh = kv(I*g] /v?)*. Una 
posible limitación en la aplicación de esta ley 
es la dependencia de d con respecto al número 
de Reynolds (1883) (Re), d > 1 cuando Re > 
0yd> 1/2 cuando Re > os, Por ejemplo, en 
el riego por gravedad es usual aportar un 
caudal constante, la presencia de un sumidero 
(infiltración del agua en el suelo) induce que el 
caudal sea nulo en el frente de onda, y puesto 
que Re «< q, es claro que Re varía desde cero 
hasta un valor máximo proporcional al caudal 
unitario de aporte. 

Existen investigaciones experimentales y 
teóricas sobre la dependencia de este exponente 
con respecto al número de Reynolds. En efecto, 
introduciendo el esfuerzo cortante en la pared 
1, = pgR,J se definen las escalas de velocidad 
(u,) y longitud (A) siguientes: 


u,= [Ty /P : h.=v/u, (1 
donde p es la densidad del agua y u. es 
denominada la velocidad de fricción. 

Con las escalas de la ecuación (1), la ley en 
potencia adquiere la forma U/u. = k(h/1J*, 
donde 0<8=2d-1<1. Ésta puede ser deducida 
aceptando que el perfil de las velocidades 
medias en potencia u(y) = k,u.(y/AJ* es válido 
en todo el dominio 0 < y < h, a partir de la 
expresión que define el caudal unitario: 


h 
q=hU¿= [u(y)dy (2) 


0 


con k, = k(1 + 0) y en donde se ha agregado el 
subíndice C, para denotar que se trata de la 
velocidad media en un canal rectangular muy 
ancho. Para el caudal en un tubo se tiene: 


R 
Q=xR"Up=Ju(y)2(R=y)dy (3) 
0 


De acuerdo con las experiencias en tubos 
llevadas a cabo por Nikuradse, los parámetros 
k, y 9 varían con el número de Reynolds 
(Re). En tubos lisos y para 5 x 10” < Re < 10”, 
Blasius propone k, = 8.74 y 8 = 1/7 (Sotelo- 
Ávila, 1974). Barenblatt (1993) y Barenblatt 
y Prostokishin (1993) argumentan que el 
comportamiento asintótico (Re —> >o) de k, y 
$ puede ser k, =k,/e + k, + k,e + ... y 6 =8,8 + 
9,e? + 8,€ + ..., donde e = 1/h(Re) y sugieren, 
a partir de los datos de Nikuradse, los valores 
empíricos k, =0.577, k, =2.50 y 9, = 1.5, y nulos 
el resto de los coeficientes. Zagarola y Smits 
(1997), con nuevos datos sugieren k, = 0.7053, 
k, = 0.3055, 6, = 1.085, 9, = 6.535, y nulos el 
resto de los coeficientes, para el rango del 
número de Reynolds de 3.1 x 10* a 3.5 x 10”. 
La dependencia de 3 con respecto al número 
Reynolds pequeños no es bien conocida. 

La aplicabilidad de la ley de resistencia en 
potencia es posiblemente limitada, ya que la 
variación del exponente en la vecindad del 
régimen laminar no es plenamente conocida. 
Ante esta dificultad se investiga la posibilidad 
de extender la ley de resistencia que resulta 
del perfil logarítmico de las velocidades de 
Prandtl a números de Reynolds pequeños. 

En el flujo turbulento cercano a la pared 
se distinguen tres regiones principales (e.g. 
Dally y Harleman, 1969; White, 1979): a) la 
región interior, donde los esfuerzos viscosos 
son dominantes, también denominada sub- 
capa viscosa; b) la región exterior, donde los 
esfuerzos turbulentos son dominantes, y Cc) 
la región intermedia o capa de solape, donde 
ambos tipos de esfuerzos son importantes. En 
la supcapa viscosa, la velocidad es una función 
básicamente lineal de la distancia a la pared: 
u/u, =y/k, cuando O < y/A, < 4. La región 
-u)/u, 
¿, es la velocidad máxima 


á: 


exterior sigue una ley de deflexión (u 
= G(y/5), donde u,, 
ubicada en el espesor de la capa límite 6, 


máx 


generalmente asimilada al radio R en un tubo. 
En la capa intermedia se ha argumentado que 
la velocidad sigue una ley logarítmica (Sotelo- 
Ávila, 1974): 
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ol (a) 
U, K |yl 


donde k y y son dos parámetros empíricos 
adimensionales. El parámetro k es conocido 
como la constante de von Kármán. 

Los valores experimentales reportados 
en la literatura se encuentran alrededor de 
los valores medios k = 0.4 y y = 0.13; algunos 
autores señala que en particular la constante de 
von Kármán puede depender del número de 
Reynolds (e.2. Daily y Harleman, 1966; Landau 
y Lifchitz, 1989). Un análisis sobre el sujeto se 
hace ulteriormente. 

La ley logarítmica es válida experimental- 
mente cuando y/A, > 30 a 70, y/8 > 0.15. En 
cuanto a la función G, en la ley de deflexión se 
ha argumentado experimentalmente que es una 


=4)/u,= —«n(y/8) 


máx 


función logarítmica: (u 
+ B, donde B es una constante. 

Considerando que la ley de deflexión puede 
ser incluida aproximadamente por la ecuación 
haciendo B = 0, y que la contribución al caudal 
por la subcapa viscosa es despreciable, la 
velocidad media se obtiene, a pesar de que la 
velocidad calculada en la pared es u(0) = -o 
en lugar de cero, como el valor principal en el 
sentido de Cauchy de las integrales que definen 
el caudal en las ecuaciones (2) y (3). La forma 
general es la siguiente: 


(5) 


donde f = ye =0.35 en el canal; PB = ye?/?= 0.58 en 
un tubo de radio R completamente lleno, h se 
reemplaza por R. 

Es claro que la ecuación (5) no es aplicable 
cuando h=0 0 R=0, ya que proporciona U(0) = 
=<o en lugar de cero. Este comportamiento de la 
velocidad media y del perfil de las velocidades 
se debe a que en la deducción de la ecuación (4) 
se desprecia la fricción debida a la viscosidad 
frente a la fricción por la turbulencia. No 
obstante, la ley de resistencia al flujo representa 


adecuadamente los datos experimentales para 
altos números de Reynolds. 

El objetivo de este trabajo es la extensión 
del perfil logarítmico de las velocidades, a 
fin de obtener el comportamiento correcto en 
la pared y su vecindad, para establecer la ley 
de resistencia, que incluya la ley de Poiseuille 
cuando el número de Reynolds sea pequeño. 


Longitud de mezcla hiperbólica 


Para alcanzar el objetivo es necesario recordar 
que la ley de Newton, válida en régimen 
laminar, establece que el esfuerzo cortante 
(1) y el gradiente vertical de las velocidades 
horizontales son proporcionales, es decir: 


(6) 


donde pv = u es el coeficiente de viscosidad 
dinámica. 

El perfil parabólico se obtiene tomando 
en cuenta que el esfuerzo cortante varía 
linealmente, siendo nulo en la superficie libre 
del agua o en el centro del tubo, y adquiriendo 
su valor máximo en la pared (1), es decir: 

T=1,(1-y/h) (7) 
el perfil de velocidades resultante es el 
siguiente: 


(8) 
Las velocidades medias en el canal y el tubo 


se obtienen introduciendo la ecuación en las 
ecuaciones (2) y (3): 


1 2 
U-=—h 9 
Cay 8) 5 
Lu 
Uy = —R 10 
T= 8 al (10) 


En régimen turbulento, el cortante, según 
las hipótesis de Reynolds, es definido por 
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la media del producto de las fluctuaciones 
alrededor de las medias de las velocidades en 
el sentido horizontal (u”) y en el sentido vertical 
(v”), es decir (Landau y Lifchitz, 1989): 


1=-p(uv”) (11) 


Para la modelación del cortante turbulento 
citemos textualmente a Levi (1989): “En la 
teoría cinética de los gases se había utilizado 
un concepto importante: el de “camino libre 
Prandtl 
una longitud característica capaz de definir el 


medio”. Análogamente, introduce 
estado de turbulencia, longitud que se podría 
interpretar como la distancia que una partícula 
fluida recorre sin perder su individualidad 
al mezclarse con sus vecinas, o bien como el 
diámetro de los agregados de partículas que 
avanzan juntas; distancia y diámetro que 
—cuando el número de Reynolds es bastante 
grande— tendrían que ser proporcionales entre 
sí, pues la energía cinética de la masa debe ser 
igual al trabajo realizado por las partículas para 
vencer la resistencia a su avance que las vecinas 
oponen. A esta longitud característica, Prandtl 
la llama longitud de mezcla” (Mischungsweg) 
y la representa con la letra /, basándose en la 
hipótesis de que, en un movimiento dirigido 
en promedio según el eje x (o sea y = 0), 
las partículas fluidas tienden a abandonar 
la capa a la cual pertenecen y —luego de 
cruzar, paralelamente al eje y, la distancia (— 
alcanzan otra posición, donde la velocidad 
media es distinta de la que la partícula poseía 
originalmente. Ahora, por ser / muy pequeña, 
la diferencia de velocidades u entre la capa 
de salida y la de llegada es / du/0y; y esta 
expresión representaría —para Prandtl— la 
fluctuación de velocidad u' = / du/dy. Como 
otras partículas se mueven de manera análoga 
en dirección opuesta, habrá choque de masas 
fluidas con valores contrarios de u'; por 
tanto nacerán fluctuaciones v” (en el sentido 
transversal) del mismo orden de magnitud 
que u'. En conclusión, el esfuerzo aparente 
-p(u'v') será proporcional a p/(du/dy)”. 
Para establecer la igualdad, Prandtl (1925) 


propone que el coeficiente de proporcionalidad 
sea absorbido por la longitud de mezcla, y 
escribir: 


du du 
dy|dy 


en la cual se ha introducido el subíndice P a la 


T=pf (12) 


longitud de mezcla de Prandtl; se ha eliminado 
la barra que denota el perfil de las velocidades 
medias; se han reemplazado las derivadas 
parciales por las totales, y se ha introducido la 
función valor absoluto para que el signo quede 
determinado por du / dy. 

El perfil logarítmico de las velocidades, 
definido por la ecuación (4), se obtiene de 
la ecuación (12), bajo el supuesto de que el 
esfuerzo cortante es una constante e igual a su 
valor en la pared (1=1, =pu.”) y que la longitud 
de mezcla en un punto dado es proporcional a la 
distancia a la pared: /, = y. Como se ha visto, la 
ecuación (12) no es aplicable en la pared debido 
a que se han despreciado los efectos viscosos. 
Lo anterior puede ser corregido, considerando 
que el esfuerzo cortante es el resultado de los 
efectos viscosos y turbulentos, es decir: 


t= Tuisc + Veurb (13) 


La sustitución de las ecuaciones (6) y (12) en 
la ecuación (13) conduce a: 


du | du 
AE (14) 


La ecuación (14) se escribe bajo la forma 
de la hipótesis de Newton, reemplazando 
el coeficiente de viscosidad cinemática en la 
ecuación (6) por un coeficiente de “viscosidad 
cinemática total' (v,), a saber (Newton, 1985): 


d 
Py (15) 


pa pee 2 
en la cual v, = v + VW... Y Van = Ldu/dy, es 
decir la “viscosidad cinemática turbulenta” es 


z , 
proporcional a u'/,. 
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De acuerdo con la ecuación (1), se define 
una longitud de mezcla total como 1 = v,/u,, de 
modo queA=1, cuando v,=v. En consecuencia, 
las ecuaciones (6), (12) y (15) se escriben de 
manera adimensional, introduciendo las 
variables adimensionales y, =y/A, u,=u/u,, T, 


=1/1,A,=2/2,=v,/v, a saber: 


du, du, 


—E o ñ == = a du, 
du * visc+ — du * turb+ P+ 
Y, Y+ 


dy, 


du. (16) 
dy, 


donde T 

A. 
Puesto que v, = vA/A, = u.A se tendrá u. 

=v+ L2du / dy, o en escritura adimensional A, 


== Viso) Ty Teurb+ 5 Taiol Ty lo, = el 


visc+ 


=1+ rd Mu, /dy,, y considerando la ecuación 
(16) se tiene 4, =1 + /?, 1, /A,. De esta última se 
deduce: 

¿=»,(A,-1) ; E = Lo,T, (17) 
donde /, = /A,, con l una longitud de mezcla 
efectiva, que coincide con la longitud de mezcla 
de Prandtl cuando 1 = T,. 

Las relaciones entre los diferentes esfuerzos 
cortantes se deducen de las ecuaciones (16) y 
(17), a saber: 


Weisc— (y %,) Te; 


Tvisc+P Veurb+ Y 


Teurbi (1-1/A,)r, , 18) 


Para 1, =1 y cuando y >> A según Prandtl 
se tiene / - ky, y de la ecuación (17) se infiere 
también A - «y (A, =1/1,>>1);en la pared (y = 
0) se tiene 2 =2, (A, =1). Se propone unir estos 
comportamientos a través de la hipérbola 
1?= 1 + by, + (xy )?. Para que existan las 
inestabilidades que originan la turbulencia se 
acepta que el perfil de velocidades presente 
un punto de inflexión a una distancia y, > 0 
de la pared (Landau y Lifchitz, 1989). 

De la ecuación (16) se deduce du /dy? = 
-(1/M dh, /dy, para 1, = 1, y en consecuencia 
dh /dy, = 0 en un punto cercano a la pared. 
La hipérbola en cuestión presenta un mínimo 


en y, = -—b/2x?, implicando b < 0. Haciendo b = 
—2xa. con o. > 0, la ecuación de la hipérbola se 
escribe como sigue: 


»,(y,)=/1-0+ (ya (19) 


El punto mínimo de la hipérbola queda 
definido por y, , =0./K, 2, = o? y se tiene 
que A, =1 tanto en la pared como a una distancia 
Y,,= 24, de la misma, por consiguiente A, <A, 


m+— 


m+ 


< 1 en el intervalo 0 < y, < ys, . 

De acuerdo con la ecuación (17), en el 
intervalo 0 < y, < Ys, , Se tiene £=-A (1-4) 
< 0, es decir que /, es imaginaria, ya que A, 
< A < 1. Esto se traduce en que el esfuerzo 
cortante turbulento definido en la ecuación 
(18) cambia de signo en esta región vecina a la 
pared. El hecho de que /, es imaginaria puede 
ser interpretado considerando que partículas 
de fluido en el plano xy, donde x es la dirección 
principal del flujo, son transferidas a otro plano 
x'y paralelo al anterior; otras partículas llegan 
al plano xy procedentes de planos paralelos, de 
modo que la masa se conserva en cada plano. 


La igualdad 1,,... = T,,,, se da cuando A, = 2 y 


corresponde a y,,, = (1+,/1+3/07 o. 


Un valor teórico de la constante de von 
Kármán 


turb+ 


Los valores de los parámetros A. y k pueden 
ser estimados bajo algunas consideraciones. 
Derivando la ecuación (17) se obtiene: 


dh, 
dy, 


(20) 


de 
£ —E=(A-—1 
a 


Ahora bien, el punto mínimo de la longi- 
tud de mezcla hiperbólica da,/ dy] =0 


Ym+ 
implica que d(. / dy | = 0. De acuerdo con la 


ecuación (20), este resultado se obtiene también 
si A,,=%%, lo que implica a1.=4/3=0.8660. De 
la ecuaciones (16) y (17) se deduce que en el 
punto crítico se tiene 2. = 


Tic. = 27, Y T 1 


turb+ + 


A O seal, = Yi, 


visc+ 
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En la región 0< y< y,, las escalas importantes 
son precisamente las que forman el punto 
crítico (y,,A,) ; la primera caracteriza una masa 
crítica y la segunda es una longitud de mezcla 
crítica que caracteriza las inestabilidades que 
originan la turbulencia. Haciendo intervenir la 
viscosidad cinemática, el análisis dimensional 
proporciona las velocidades correspondientes 
u, =VW/y, y v, =v/A,, así como el cortante en la 
pared T, = Pu, Di este cortante tiene la forma de 
aquel definido por la ecuación (11). Dado que 
T, = pu”, ecuación (D), se tiene u?=u,0,, y M = 
y,,A,, es decir, en escritura adimensional u,, 0... 
=1 y y, A, = 1. De modo que A, = Y implica 
que y, = 0/k =2, esto es k=4/3=0.4330 y, 
ver ecuación (23), y=2/3-1=0.1547, valores 
muy cercanos a los reportados en la ecuación 
(4). Por la importancia de estos parámetros se 
establecen a continuación los valores resultan- 
tes del análisis precedente: 


0.=3/3= 0.8660;  k=v3=0.4330; 
(21) 


y =3v3 -1=0.1547 


Es importante señalar que este valor teórico 
de la constante de von Kármán es muy cercano 
al valor obtenido por Zaragola (1996), a partir 
de ensayos realizados en el rango 3.1 x 10*a 3.5 
x 10” del número de Reynolds más amplio que 
el rango estudiado por Nikuradse (1932, 1933) 
(Zaragola y Smits, 1997; McKeon et al., 2004), a 
saber: k= 0.436. 

Con los valores de ar y k definidos por la 
ecuación (21) se deduce la frontera superior de 
la región crítica cercana a la pared y, = 4 y la 
distancia a la pared donde los esfuerzos viscoso 
=2(1+45)= 6.4721. 
Es de notar que el valor de y,, = 4 coincide con 


y turbulento son iguales y,,,. 
el valor experimental del espesor de la subcapa 
viscosa. 

Las gráficas de la longitud de mezcla 
hiperbólica definida por la ecuación (19) 
y la asíntota «y, — Q, con los valores de los 
parámetros proporcionados por la ecuación 
(21), se muestran en la figura 1. 


Distancia a la pared 


0 1 2 3 4 
Longitud de mezcla 


Figura 1. La longitud de mezcla hiperbólica, ecuación 
(19), y su asíntota. Los valores de las constantes son los 
definidos por la ecuación (21). 


Perfil de las velocidades 


El perfil de 
correspondiente a T,= 1 se obtiene de las 
ecuaciones (16) y (19), con la condición u (0) = 
0, a saber: 


velocidades  adimensional 


1, [| xy,- 
100 no] . 


en donde se ha introducido la notación u_, (y) 
para indicar que T, = 1 se obtiene de la ecuación 
(8) cuando h, > oo, 

Cuando y, >>, O la ecuación (22) tiene el 
comportamiento de la ecuación (4): 


Jl 1-0 
0.0 la de |, 2 a) 


y cuando y, > 0, su desarrollo es el siguiente: 


uo (y)= 9, +3 ox +0 (y7) (04) 
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el cual muestra el comportamiento lineal en 
la subcapa viscosa cerca de la pared (Panton, 
1984). 

El perfil adimensional de las velocidades 
correspondiente a la ecuación (7), 1,=1-y,/h,, 
se obtiene de las ecuaciones (16) y (19), con la 
condición u (0) =0, a saber: 


1 104 Ky,— 
AA E Jn : 
: 


2 Pd] 


Es importante señalar que cuando y, es 
muy grande y en consecuencia también h,, ya 
que 0 < y, < h,, de la ecuación (25) se deduce 
(1 - a/xh) 
In(y,/y). Escribiendo la ecuación (4) como 
u, = kIn(y, /y), en donde se ha cambiado la 
notación para el parámetro original de von 


el comportamiento u (y ,h,) = 


Kármán (x,,); es claro que k;! = k" (1 — a/ kh), 
la cual puede justificar la dependencia de dicho 
parámetro con respecto al número de Reynolds 
definido por Re” =h,=u.h/v, argumentada por 
algunos autores cuando el perfil logarítmico 
es aplicado en flujos donde el esfuerzo 
cortante no puede ser considerado constante. 
Claramente, cuando h, > es, la ecuación (25) se 
reduce a la ecuación (22). 

El comportamiento de u, (y ,h,) cuando y, > 
0 es el siguiente: 


2 
u,(yyh,)=y,-(- o +0 (y?) (26) 
+ 


en donde de nuevo aparece el comportamiento 
lineal en la subcapa viscosa. 

El perfil parabólico correspondiente al 
régimen laminar y proporcionado por la 
ecuación (8) en escritura adimensional es u, = 
y, - y? /2h_,el cual se deduce de la ecuación 
(26) cuando h,>0. 

La velocidad máxima ocurre a la distancia 
y, = h, de la pared; la diferencia entre esta 
velocidad y aquella en y, cercana ah, se calcula 
con las ecuaciones (22) y (25). Los resultados 
son las leyes de deflexión siguientes: 


(27) 


1, (h, 
Ue (,)= Uso (y,) dé a) 


ES 


1 
U, (1,,h,) =U, (Y, h,) dd dé ' 


1 
ERE y,) 


donde la primera se deduce de la segunda 
haciendo h, >> y h _-y,>0. 

Los perfiles de las velocidades definidos por 
la ecuación (22) y (25) se muestran en la figura 
2, en los cuales se debe observar la existencia 
del punto de inflexión. 


Velocidades medias 


Se obtienen a continuación las expresiones 
de la velocidad media del agua en un tubo 
completamente lleno y en un canal rectangular 
muy ancho. 


Tubo. La velocidad media adimensional para T, 
= 1 se obtiene introduciendo la ecuación (22) en 
la ecuación (3): 


Distancia a la pared 


0 2 4 6 8 10 
Velocidad 


Figura 2. El perfil de las velocidades, ecuaciones (22) 
y (25), considerando esfuerzos cortantes constante y lineal. 
Los valores de las constantes son los definidos por la 
ecuación (21). 
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0? R,-0a+A, (R 
ps) eat 
KR -0 (29) 
—[(4R, +30), (R,)-3 
erario Sid 0.) Gl aja o] 
2 
=> A, (R,)-1 
ERP R)-1] 
Los comportamientos de U, (R,), cuando 
R,>0 y R,> o son, respectivamente, los 
siguientes: 


1 
Unos (R)=GR.+O(R:) (30) 
R, 


1 
Ur... (Ry) a Ea 


la ecuación (31) corresponde a la ecuación (5). 
La velocidad media adimensional para 1. = 
1 - y./R. se deduce de las ecuaciones (3) y (25): 


ll 104 (04 
R,)==|1 211 
Ur le Sl al | En) 
o ne] 
2 RE 1-0 


(A EA 
DER KR, KR 
(32) 


+ 


(31) 


OS 


ol ol ol 
| A a 
El ES, el a | 


+ 


Los comportamientos de U, (R) cuando 
R,>0 y R,> o son, respectivamente, los 
siguientes: 


Uy, (R,)=3R,+0(Ri) (33) 


1 R 
Ur, (R, = ed (os ) (34) 


En escritura adimensional, la ley de 
Poiseuille de los tubos capilares definida por la 
ecuación (10) es U, (R ,) = UR, la cual se deduce 


precisamente de la ecuación (33). 


velocidad media 
adimensional para Tt, = 1 se obtiene introdu- 


Canal rectangular. La 


ciendo la ecuación (22) en la ecuación (2): 


Uc(h)= ad [e 0) 


«h 1-0 


(85) 
1 
E [a (1,)-1] 


Los comportamientos de U,, (hn) cuando 
h,=> 0 y h, > « son, respectivamente, los 
siguientes: 


Uc., (t1,)=3h, +0 (17) (36) 
A 
Uco (1,) ==In ( (37) 


la ecuación (37) corresponde a la ecuación (5). 

La velocidad media adimensional para 1, = 
1 — y, /h, se obtiene introduciendo la ecuación 
(25) en la ecuación (2): 


1-a 


(53) p Pato 


1-02 ES a+), eo] (68) 


Los comportamientos de U, (h,) cuando 
h,=> 0 y h, > « son, respectivamente, los 
siguientes: 
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Ue, (1,)=3h,+0 (12) (39) 


1 h 
Uc, (,) 5 o () (40) 


En escritura dimensional, la ecuación (39) 
corresponde al régimen de Poiseuille sobre una 
placa plana, ecuación (9). 


El factor de fricción 


La teoría anterior es válida en paredes lisas. 
En paredes rugosas y a grandes números de 
Reynolds, la escala A, se reemplaza por la escala 
A, que caracteriza la rugosidad de la pared; 
siguiendo a Colebrook y White (Colebrook, 
1939), en régimen intermedio, la escala A. Se 
reemplaza por la escala A, definida como A, = 
A, +A, 

A altos números de Reynolds, la escala 
A, es proporcional a la altura media de las 
rugosidades de la pared, tomada como la altura 
equivalente de granos de arena (e) o escala de 
Nikuradse, es decir: 


Jim A, (Re, £)=c8 (41) 


donde c es un coeficiente de proporcionalidad. 

En un tubo, la ley de resistencia se aplica 
utilizando el coeficiente de fricción de Fanning 
(f) definido en la ecuación de Darcy-Weisbach 
J = fu? / 22D, donde D es el diámetro interno 
del tubo (Darcy, 1857; Darcy y Bazin, 1865; 
Weisbach, 1850). Considerando la velocidad 
de fricción, este coeficiente queda definido 
como f = 8/U,?. Reemplazando la escala A, 
por la escala A, = A, + A, la relación entre el 
coeficiente de fricción, el número de Reynolds 
definido por Re = UD/y, y la escala relativa 
definida por A,, = 2,/ D, se obtiene de la 
ecuación (29) o (32): 


O 1] e 
nllatyal] 


donde ),, = ce, cuando Re > co y e, =8/D es la 
rugosidad relativa de Nikuradse. 

El parámetro c se estima a partir de la 
información obtenida por Nikuradse en sus 
experiencias realizadas en tubos con seis 
valores de rugosidad relativa a grandes 
números de Reynolds (Sotelo-Ávila, 1974). 
Los valores de (e/D,f) se presentan en el 
cuadro 1; asimismo, se reportan los valores 
correspondientes de A, / D obtenidos utilizando 
en la ecuación (42) las ecuaciones (29) y (32) 
con Re > oo. La relación entre 2,/D y e/D se 
muestra en la figura 3. Para la ecuación (29), c 
= 0.125 con R? = 0.9996, mientras que para la 
ecuación (32), c = 0.09 con R? = 0.9997. 

En la zona de transición, la función A. (Re, 
e) es desconocida y podrá ser aproximada 
empíricamente a partir de los datos experi- 
mentales. Con la finalidad de que la función 
A[Re, 1. (Re, e)l, que define al factor de fricción 
implícitamente a través de la ecuación (42), 
manifieste algunas de sus características, se 
tomará el comportamiento asintótico definido 
por la ecuación (41). 

La función f(Re,ce), mostrada en la figura 4 
para los seis valores de la rugosidad relativa 
de Nikuradse, utiliza la ecuación (32). La 
función así construida no proporciona el com- 
portamiento observado del factor de fricción 
en la zona de transición; su comportamiento 
asintótico es similar al del factor de fricción 
mostrado en el clásico diagrama de Moody 
(1944, 1947). Sin embargo, el diagrama de 
la figura 4 tiene la principal ventaja de que 
contiene como caso límite, a bajos números 
de Reynolds, al régimen laminar de Poiseuille 
para el cual f=64/Re. 

En el caso de tubos o canales rugosos, si la 
rugosidad de la pared se obtiene de manera 
inversa en un problema de flujo en régimen 
laminar o de transición, es posible que el 
valor no tenga ningún significado físico, por 
ejemplo, que sea del mismo orden que el del 
diámetro del tubo, debido a que f(Re, ce) es 
aproximada en esta región de flujo. En esta 
situación, tal vez es mejor retener el valor 
obtenido de la escala 1 (Re, e), la cual estará 
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Figura 3. Relación entre la escala de rugosidad en la pared (A.,) relativa al diámetro del tubo (D), y la rugosidad 
relativa de los experimentos de Nikuradse (e/D), ecuación (42): a) con la ecuación (29) (A/D) = 0.125(€/ D), 
con R?= 0.9996 ; b) con la ecuación (32) (A/D) = 0.09(€/D), con R?= 0.9997. 


Cuadro 1. La rugosidad relativa (e/D) y el factor de fricción de Fanning (f), correspondientes a las experiencias 
de Nikuradse a grandes números de Reynolds (Re > 09); los valores de fson obtenidos de las gráficas de 
Sotelo-Ávila (1974). La escala relativa A, /D es obtenida a partir de las ecuaciones (29) y (32) en la ecuación (42). 


(A/D) x 10* (A/D) x10* 
e/D Í ecuaciones (29) ecuaciones 
y (42) (32) y 42) 
1 
30 0.0600 41.84 29.85 
1 
612 0.0450 20.58 14.70 
1 
1720 0.0350 9.88 7.06 
1 
252 0.0275 4.37 SiS) 
1 
504 0.0235 2.41 1.73 
1 
1014 0.0200 1.24 0.89 


asociada con un número de Reynolds del 
régimen de flujo y con una rugosidad de la 
pared, y ser utilizada como un valor medio 
en un pequeño rango de variación de este 
número. 


Conclusiones 


Se ha considerado que el esfuerzo cortante en 
el flujo del agua es el resultado de la suma 


de un esfuerzo viscoso y de un esfuerzo 
turbulento. El primero es proporcional, de 
acuerdo con Newton, a la derivada del perfil 
de las velocidades, mientras que el segundo, 
de acuerdo con Prandtl, al cuadrado del 
producto de la derivada de este perfil y una 
longitud de mezcla. Se ha introducido una 
longitud de mezcla hiperbólica para escribir 
la relación entre el esfuerzo cortante total y la 
derivada del perfil en la forma de Newton. 
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Figura 4. El coeficiente de fricción de Fanning (f) en un tubo completamente lleno de agua, ecuación (42) 
considerando la ecuación (32), en función del número de Reynolds (Re) en un tubo liso y para los seis valores 
de la rugosidad relativa de Nikuradse (€ /D). 


La longitud hiperbólica ha sido construida 
de tal manera que tenga el comportamiento 
lineal con respecto a la distancia a la pared 
de la de Prandtl, en la cual el coeficiente de 
proporcionalidad es la constante de von 
Kármán, y que tenga un valor constante en 
la pared para incluir el régimen laminar. Esto 
induce que la hipérbola presente un punto 
mínimo cerca de la pared, hecho que a su vez 
induce un punto de inflexión en el perfil de 
las velocidades. La existencia de este punto 
crítico permite deducir un valor teórico de la 
constante de von Kármán bastante cercano a 
su valor experimental. 

Con el perfil de las velocidades obtenido 
se han deducido las fórmulas para calcular 
las velocidades medias del fluido en un 
tubo completamente lleno de agua y en un 
canal rectangular de ancho muy grande; 
a partir de estas fórmulas se obtienen las 
leyes de resistencia al flujo. Se ha estudiado 
en particular la ley de resistencia en el tubo 


llevándola a la forma de la ecuación de Darcy- 
Weisbach; el factor de fricción resultante 
presenta los comportamientos extremos 
clásicos, en función del número de Reynolds, 
es decir el régimen laminar de Poiseuille a 
bajos valores y el de Prandtl a altos valores de 
este número. 

El efecto de la rugosidad de la pared a 
altos números de Reynolds se ha introducido 
en la ley de resistencia, siguiendo la idea de 
Colebrook y White de sumar a la escala de 
la turbulencia una escala proporcional al 
tamaño de los granos de arena de Nikuradse. 
El factor de fricción resultante tiene el mismo 
comportamiento que el mostrado en el 
diagrama clásico de Moody. Sin embargo, 
esta escala depende del número de Reynolds 
en régimen transitorio de acuerdo con los 
resultados de Nikuradse. 
La dependencia de esta escala a pequeños y 


experimentales 


moderados valores del número de Reynolds y 
de la rugosidad es un sujeto de investigación. 
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Abstract 


FUENTES, C. Hyperbolic mixture length, velocity profile, and flow resistance law. Water 
Technology and Sciences, formerly Hydraulic engineering in Mexico (in Spanish). Vol. 
II, No. 1, January-March, 2011, pp. 19-31. 


Considering that in water flow total shear stress is the sum of a viscous shear stress and a 
turbulent shear stress, a hyperbolic mixture length is introduced in the Newton hypothesis 
to describe the relationship between total shear stress and the mean velocity profile derivative. 
The hyperbolic length is reduced to the Prandtl length which is proportional to the wall 
distance when the latter is large. The existence of an inflexion point in the velocity profile 
induces a theoretical value of the von Kármán constant very close to its experimental value. 
From this, profile formulas for average fluid velocities in a pipe and rectangular channel are 
deduced, and from these latter, resistance to flow laws are established. The pipe resistance 
law, expressed in a Darcy-Weisbach form, indicates that the friction factor shows the extreme 
behaviors of Poiseuille and Prandtl laws. The wall roughness for high Reynolds numbers is 
introduced following the Colebrook and White idea, adding to the turbulence scale another one 
proportional to the Nikuradse sand grains size; the friction factor has the same behavior as that 
of the Moody diagram. The last scale will depend on the Reynolds number in transient regime 
in agreement with the experimental results of Nikuradse and its dependency with respect to 
small and moderate Reynolds number values, and the wall roughness is a study subject. 


Keywords: theoretical value of the von Kármán constant, wall roughness effects. 
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